Anisotropy in Born-Infeld brane cosmology 
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The accelerated expansion of the universe together with its present day isotropy has posed an 
interesting challenge to the numerous model theories presented over the years to describe them. In 
this paper, we address the above questions in the context of a brane-world model where the universe 
is filled with a Born-Infeld matter. We show that in such a model, the universe evolves from a highly 
anisotropic state to its present isotropic form which has entered an accelerated expanding phase. 
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I. INTRODUCTION 



the FRW model. 
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The problem of the diverging Coulomb field and self- 
energy of point particles in Maxwell's theory of electrody- 
namics led Born and Infold (BI) in 1934 to the construc- 
tion of a nonlinear extension of classical electrodynam- 
ics [1]. In recent years, the study and use of nonlinear 
electrodynamics has been the focus of much attention in 
the context of cosmology. Such strong motivations stem 
from developments in string/M-theory, where it has been 
shown that the BI theory naturally arises in the low en- 
ergy limit of the open string theory [2]. Nonlinear elec- 
trodynamics has also been applied to several branches of 
physics, namely, the effective theories at different levels 
of string/M-theory [3], cosmological models [4, 5], black 
holes [6, 7] and wormhole physics [8], among others. 

The ubiquitous Friedmann-Robertson- Walker (FRW) 
model has long been considered standard against 
which other model theories are compared. In this model 
the universe is regarded as homogeneous and isotropic 
with only one dynamical feature, its rate of expansion 
or contraction. Its high degree of symmetry makes it 
easy to analyse but also makes it somewhat unrealistic. 
However, it does model the average properties of the ob- 
served universe, on the largest scales, quite well. The 
Bianchi models, on the other hand, pose a more realis- 
tic scenario. These models represent a homogeneous but 
anisotropic universe which are essentially indistinguish- 
able from FRW models far from the initial singularity 
or, in the case of a collapsing universe, far from the fu- 
ture singularity. The present epoch is far from all such 
singularities. Bianchi models have three dynamical pa- 
rameters, the expansion or contraction in the three spa- 
tial directions which are, in general, all different. The 
three rates of expansion or contraction asymptotically 
(in time) equalize, thus producing a model which looks 
like that of a FRW. If the universe did not begin by ex- 
panding exactly isotropically, then the Bianchi models 
would describe its early behavior more accurately than 
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The idea that our four-dimensional (4Z?) universe is a 
hypersurface (brane) in a 5D spacetimc (bulk) [9-11] has 
been the dominant idea behind many model theories at- 
tempting to explain the observed universe over the past 
decade. One of the most successful of such higher dimen- 
sional models is that proposed by Randall and Sundrum 
where the bulk has the geometry of an AdS space admit- 
ting Z2 symmetry [10]. They were successful in explain- 
ing what is known as the hierarchy problem; the enor- 
mous disparity between the strength of the fundamental 
forces. The Randall-Sundrum (RS) scenario has had a 
great impact on our understanding of the universe and 
has brought higher dimensional gravitational theories to 
the fore. In certain RS type models, all matter and gauge 
interactions reside on the brane while gravity can prop- 
agate into the bulk. In obtaining the field equations on 
the brane one uses the Israel junction conditions [12] and 
the Gauss-Codazzi equations, as employed by Shiromizu, 
Maeda and Sasaki (SMS) [13]. These field equations dif- 
fer from the standard Einstein field equations in 4D in 
that they have additional terms like tt^ v which depend 
on the energy-momentum tensor on the brane and the 
electric part of the Weyl tensor E^ v , leading to the ap- 
pearance of a quadratic term in the Friedmann equation. 
This term was initially considered as a possible solution 
to the accelerated expansion of the universe. However, 
soon it was realized to be incompatible with the big bang 
nucleosynthesis, requiring additional fixes [14]. 

Since their introduction, the RS type models have been 
employed frequently in conjunction with a myriad of mat- 
ter fields to study the behavior of the observed universe. 
What has been lacking however is such a study with non- 
linear electrodynamics of the BI as the matter field. How- 
ever, as we will show in the Appendix, the BI matter can- 
not be used in a brane-world scenario together with Is- 
rael junction conditions. The reason is that the solutions 
of the resulting field equations predict a universe which 
is incompatible with present day observations. In addi- 
tion, there have been arguments on the uniqueness of the 
junction conditions [15] or their use when more than one 
non-compact extra dimension is involved. In view of the 
above, one is lead to consider an alternative brane-world 
scenario in which the Israel junction conditions become 
redundant. Brane-world scenarios under more general 
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conditions and still compatible with the branc- world pro- 
gram have therefore been rather extensively studied over 
the recent past where it has been shown that it is possi- 
ble to find a set of cosmological solutions in accordance 
with the current observations. Under these conditions, 
without using Z 2 symmetry or postulating any junction 
conditions, the Friedmann equation is modified by a ge- 
ometrical term which is defined in terms of the extrinsic 
curvature, leading to a geometrical interpretation of dark 
energy [16] . An example of such a scenario was presented 
in [17] where particles are trapped on a AD hypcrsurface 
by the action of a confining potential. The dynamics of 
test particles confined to a brane by the action of such a 
potential at the classical and quantum levels were stud- 
ied in [18]. In [19] the authors showed that the quadratic 
term in the energy-density can be produced in this model, 
suggesting that the model is identical to the SMS method 
at least in the isotropic universe. In this paper, we con- 
sider an anisotropic brane-world with Bianchi type I ge- 
ometry filled with a nonlinear electromagnetic field of the 
BI type. The electromagnetic tensor has six independent 
components of which only two survive to construct the 
BI energy-momentum tensor. The anisotropy parameter 
and the deceleration parameter are calculated and shown 
to be compatible with present observations, namely that 
the universe is in an accelerating phase and expanding 
isotropically. 



scenario a new conserved quantity, Q^ u , appears which is 
totally geometric in nature and depends on the extrinsic 
curvature only, reflecting the effects of the extra dimen- 
sion. In the SMS method one replaces this tensor with 
a tensor build up by the energy-momentum tensor using 
the Israel junction conditions. We will compute this ten- 
sor using the Codazzi equation. The field equation of the 
model is given by [20] 
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where a is the energy scale of the brane, A is the cos- 
mological constant, is the electric part of the Weyl 
tensor and Q M „ is a conserved quantity, V \iQ ilv = 0, ex- 
pressed in terms of the extrinsic curvature and its 
trace 

Q MV = KK MV - K^K" + i (K a(s K afi - K 2 ) 9fi „. (5) 
In order to calculate we first use the York's relation 
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where £ is a vector normal to the brane. This tells us that 
the extrinsic curvature is diagonal and does not depend 
on the spatial coordinates. After separating the spatial 
components, the Codazzi equation relates the diagonal 
elements of the extrinsic curvature tensor [201 



II. THE SETUP 

Let us start by defining the Bianchi type I universe 
through the metric 



ds 2 = -dt ®dt + ^2 aiitfdx 1 - (g> dx l 



(1) 



where Oj are the scale factors. We also have the usual 
definition of the volume scale factor v, the directional 
Hubble parameters Hi and the mean Hubble parameter 
H, defined as 

3 . 3 

u = TTa i , Fi = -, m = Y^H h 

l—l 1=1 

AH i =H i -H. 1 = 1,2,3. (2) 

The physical observables are the anisotropy parameter A 
and the deceleration parameter q given by 



, = - ff --l. (3) 



^ j,k ~ ^ % k,j-i hjjk — 1)2,3, 

K\ + —K\ = —K° . (no summation) (7) 
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The first equation simply emphasizes the result that the 
spatial components of the extrinsic curvature do not de- 
pend on the spatial coordinates. The choice of the ansatz 
A'oo = d(t) for the temporal component of the extrinsic 
curvature and use of the second equation in (7) leads one 
to the other components of the extrinsic curvature 
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with Cj's being constants. Using equation (5) we can 
obtain the components of the tensor Q^ v 

Q°o = - (eie 2 + e 2 e 3 + e 3 ei) 

( c l c 2 c 3 / 
— (e2 + e 3 ) H [ex + e 3 ) H (e 2 + e 3 ) 
ax a 2 a 3 

ax a 2 a 2 a 3 a 3 ax 
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In view of the discussion above, we focus attention on 
a brane world scenario [18] in which the use of Z 2 sym- 
metry and Israel junction conditions are relaxed in favor 
of a confining potential V whose role is to localize the 
gauge fields of the standard model on the brane. In this 
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The components Q 2 2 and Q 3 3 are obtained from (11) by 
permutation of indices. Let us now concentrate on the 
case where the temporal component of the extrinsic cur- 
vature is constant, taking it as d(t) = 1. Such a choice 
can be understood in terms of the definition of and 
observation that goo = — 1. In this case the components 
of Q^ v take the simpler form 



We note that in the weak field limit P^ v s 
energy-momentum tensor can be written as 
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and again Q 2 2 and Q 3 3 are obtained by permutating the 
indices. If we assume that our 51? bulk has a constant 
curvature, £^ v vanishes and the field equation becomes 



(13) 



Let us now assume that the AD universe is filled with 
the BI field and is described by the following action 



S F = J d 4 xC(F), 



(14) 



where 

C(F) = V~gL{F) 



= D { V=5 - yj\ det ( 9fiu + p^F^) \ \ . (15) 

Here, /3 is a coupling constant with the dimension of the 
field. In four space-time dimensions equation (15) can be 
expanded to give [21] 



2 




16/3 



(16) 



(17) 



In general, the anti-symmetric tensor F^ has six inde- 
pendent components. However, since the Einstein and 
Qp, v tensors are diagonal, the six off-diagonal components 
of (13) result in six algebraic equations with solutions 
where only two non-zero components, F tx and F yz , sur- 
vive. Using the field equation for the BI matter (18) 
and conservation equation for the BI energy-momentum 
tensor, V^t^ v — 0, one can find the two remaining com- 
ponents 
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wherc fi are constants and 
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This relation ensures that the conservation equation and 
the field equation of the BI matter are satisfied and elim- 
inates the constants from the energy-momentum tensor 
which can be written using equation (20) as 
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The energy-momentum tensor of the BI field can 
be written in the form of a perfect fluid t£ = 
diag(— /0,p" ,p ± ,p ± ) with the equation of state of the form 
pll = — p. Another equation of state relating the trans- 
verse pressure can be obtained from (23) 



In this equation, F^ v = A^ v — A v>fl is the usual elec- 
tromagnetic tensor with 4-potcntial A 1 * and F^ = 
ie^^Fpt,. This action reduces to that of the Maxwell 
in the limit j3 — > oo. The field equation for the BI field 
is obtained from the above Lagrangian 



0, 



(18) 



P = 



P 



l - Ho' 



(25) 



In the limit j3 — > oo one can see from equation (25) that 
the energy-momentum tensor of the BI field takes the 
form diag(— p, — p, p, p). As can be expected, this is the 
energy-momentum of the Maxwell's field with 



where we have introduced the second rank tensor 



pllv = 1 dL(F) 
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where v is defined in (2). Equations (12), (23) and the 
Einstein tensor for the Bianchi I metric suggest that the 
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scale factors in the (yy) and (zz) directions are propor- 
tional. Solving the (yy) and (zz) components of the field 
equation (13) one obtains 

asjt) = a 2 (t) , 27 n 

C3 C 2 



This means that the second and third components of the 
Einstein equation are identical. Equation (27) also im- 
plies that the constants c 2 and C3 must have the same 
sign. It only remains to solve the field equations (13) 
which, using equation (2) can be written as 
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where A takes the form 
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Equation (29) can be integrated to find the scale factor a 2 implicitly 

da 2 
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One can also find a relation between the scale factors 
a\(t) and a 2 (t) 
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(33) 



where C5 is an integration constant. 



using equation (33). These equations show that at late 
times the BI field differ from that of the Maxwell in the 
second terms which is of order a^ 8 . Consequently, the 
BI field very rapidly turns into the Maxwell field unless 

the constant -^S^ is very large. One can also see that the 

anisotropy and deceleration parameters behave according 
to 



III. LATE TIME BEHAVIOR 

An explicit form for the solution of the above equa- 
tions seems to be unavailable. It is therefore desirable to 
analyze the asymptotic behavior of the observable quan- 
tities at late times. Assuming (3 < 00, the behavior of 
the energy-density and pressure of the BI field is found 
to be 
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Equations (34) are obtained by writing H 2 in terms of a 2 
from equation (32) followed by writing H\ in terms of a 2 
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The asymptotic behavior of the scale factor is 



a x ~ -|V3A + 9a 2 . (38) 

One can see from equations (27) and (38) that the Hub- 
ble parameter of the three scale factors become equal 
at late times, predicting an isotropic universe. In fig- 
ures (1), (2) and (3) we have plotted the anisotropy 
and deceleration parameters for three different values 
of (3, including the Maxwell's regime /3 = oo and con- 
stants c\ = —100, C2 = —1, C3 = —60. These fig- 
ures show that the behavior of the anisotropy parameter 
and the deceleration parameter are largely dependent on 
the Born-Infeld coupling. As the constant /3 decreases, 
the early time anisotropy of the BI universe becomes 
less pronounced. From the first figure we see that the 
slope of the anisotropy parameter decreases, leading to a 
less anisotropic universe which becomes isotropic with a 
smaller rate than the Maxwell universe. For (3 = 10~ 5 , 
the universe reaches an isotropic stage at early times, 
becoming less anisotropic again and finally reaches the 
present isotropic state. The deceleration parameter dia- 
gram shows the same behavior; the BI universe starts 
from a smaller amount of deceleration at early times 
and becomes accelerating considerably slower than the 
Maxwell case at late times. 

To answer the question of what value of (3 is the most 
suitable, one requires comparison against the observa- 
tional data which may become a realistic possibility in 
not too distant a future. In figure (2) we have plotted 
the early time behavior of the deceleration parameter. 
One can see from the figure that the BI universe with 
the coupling f3 = 10 -5 has an early accelerating stage 
which alternately becomes decelerating and accelerating 
at late times. Finally we note that the late time behavior 
of the theory discussed in this section does not depend on 
the values of the constants of the model. Indeed at late 
times the model asymptotes the Maxwell's theory, which 
can be seen from the figures. However the detailed early 
time behavior of the model depends on the values of the 
constants of the model, but its qualitative behavior is the 
same as the constants assume different values. 



IV. CONCLUSIONS AND REMARKS 

The observation of the present acceleration of the 
universe is one of the most intriguing problems in the 
present-day cosmology. Explanations based on the exis- 
tence of dark energy or modified theories of gravity are 
common approaches to describe it. 

In this paper we have dealt with this problem in the 
context of a brane- world model in which the assumption 
of the Z2 symmetry and use of the Israel junction condi- 
tions are relaxed and the matter content of the universe is 
assumed to be of the BI type. To calculate the conserved 
quantity appearing in the field equations, one needs 
to calculate the components of the extrinsic curvature. 
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FIG. 1: The Anisotropy parameter in the cases: /3 = 00 
(dashed-dotted line), /3 — 1 (dashed line), /3 = 1CF 5 (solid 
line) 
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FIG. 2: The early time behavior of the deceleration parameter 
in the cases: p — 00 (dashed-dotted line), (3 = 1 (dashed line), 
P = 1(T 5 (solid) 
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FIG. 3: The late time behavior of the deceleration parameter 
in the cases: /3 — 00 (dashed-dotted line), f3 = 1 (dashed 
line), = 10~ 5 (solid) 
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We have used the Codazzi equation to achieve this. The 
components of the extrinsic curvature are obtained by 
the choice of an ansatz for its (00) component. However, 
as long as the bulk metric is not specified, this component 
can be chosen in such a way as to make the tensor Q^ v 
taking as simple a form as possible. In the SMS method, 
this is achieved by the form of the energy-momentum 
tensor. However, as is discussed in the Appendix, the 
form of the BI energy-momentum tensor is not arbitrary 
and is calculated from the BI field equations, indicat- 
ing that the SMS procedure may not provide solutions 
compatible with the present observations when confin- 
ing the BI field as the matter source to the brane. In 
contrast, the method used in this paper provides more 
freedom in choosing the extrinsic curvature, making the 
confinement of the BI matter on the brane a possibil- 
ity. In other words, in the SMS method one restricts the 
bulk geometry to be satisfied by the Israel junction con- 
ditions, leaving the energy-momentum tensor arbitrary. 
In our model however, the procedure does not have such 
a restriction. 

The asymptotic behavior of the BI matter was also 
studied. As can be seen from equation (34), its behavior 
mimics that of the Maxwell field in a certain limit; the BI 
field makes the transition of the universe from a highly 
anisotropic early state to its present isotropic form slower 
than the Maxwell field and this depends on the parameter 
/3, as is seen from figure (1). 

V. APPENDIX 

Let us now discuss the consequences of using the BI 
field in brane world scenarios where the confinement of 



gauge fields is achieved through the SMS method. As 
is well known, in this method one uses the Israel junc- 
tion conditions to replace the extrinsic curvature tensor 
with the energy-momentum tensor. It follows that a new 
tensor, n^" , emerges which is quadratic in the energy- 
momentum tensor. This tensor can be written as 

(39) 

where r = r a a and g^ v is the brane metric. One of the 
consequences of the SMS field equations is that for the 
constant curvature bulk which we have considered in this 
paper, the above tensor must be conserved V^tt^ = 
[13]. If we consider the Bianchi I metric and note that 
the energy-momentum tensor for this space-time takes 
the form = diag(— p, —p,p,p), equation (23), the con- 
servation equation for -k^ v reduces to 

(p + p)(p + 2p) =0. (40) 

The first term cannot be zero as one can see easily from 
definitions of p = Tq and p = rf from equations (23). So 
the only equation of state for the BI matter, consistent 
with the above constraint is p+2p = const. This equation 
can be solved for the scale factors, with solution 

a 2 oc — . (41) 

This solution is obviously ruled out by present observa- 
tions. 
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